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To understand the nematicity in Fe-based superconductors, nontrivial ^-dependence of the or¬ 
bital polarization {AExz{k), AEyz{k)) in the nematic phase, such as the sign reversal of the orbital 
splitting between F- and X,Y-points in FeSe, provides significant information. To solve this prob¬ 
lem, we study the spontaneous symmetry breaking with respect to the orbital polarization and 
spin susceptibility self-consistently. In FeSe, due to the sign-reversing orbital order, the hole- and 
electron-pockets are elongated along the ky- and /ca^-axes respectively, consistently with experiments. 
In addition, an electron-pocket splits into two Dirac cone Fermi pockets with increasing the orbital 
polarization. The orbital-order in Fe-based superconductors originates from the strong positive 
feedback between the nematic orbital order and spin susceptibility. 

PACS numbers: 74.70.Xa, 75.25.Dk, 74.20.Pq 


The spontaneous symmetry breaking from C 4 - to C 2 - 
symmetry, so called the electronic nematic transition, is 
one of the fundamental unsolved electronic properties in 
Fe-based superconductors. To explain this nematicity, 
both the spin-nematic scenario [1-6] and the orbital order 
scenario [7-13] have been studied intensively. Above the 
structural transition temperatures Tgtr, large enhance¬ 
ment of the electronic nematic susceptibility predicted 
by both scenarios [1, 10] is actually observed by the 
measurements of the softening of the shear modulus Cqq 
[1, 10, 14, 15], Raman spectroscopy, [16-18], and in-plane 
resistivity anisotropy Ap [19]. 

To investigate the origin of the nematicity, FeSe (Tc = 
9 K) is a favorable system since the electronic nematic 
state without magnetization is realized below Tgtr = 90 
K down to 0 K. Above Tgtr, the antiferromagnetic fluc¬ 
tuations is weak and T-independent according to the 
NMR [20, 21] and neutron scattering [22-24] studies, 
in contrast to the sizable spin fluctuations above Tgtr in 
LaFeAsO [25] and BaFe 2 As 2 [26]. This fact means that 
the magnetic instability is not a necessary condition for 
the electronic nematic state. In contrast to the smallness 
of the spin fluctuations, large nematic susceptibility is 
measured by Cqq and Ap in FeSe. Based on the orbital- 
spin fluctuation theory, called the self-consistent vertex- 
correction (SC-VC) theory, the development of the strong 
orbital fluctuations in FeSe are explained even when the 
spin-fluctuations are very small, consistently with exper¬ 
imental reports in FeSe [27]: The strong orbital fluctu¬ 
ations originate from the Aslamazov-Larkin vertex cor¬ 
rection (AL-VC) that describes the orbital-spin mode- 
coupling [10]. The nematic CDW in cuprate supercon¬ 
ductor also originates from the AL-VC [28, 29]. 

The nontrivial electronic state below Tgtr gives a crucial 
test for the theories proposed so far. In the orthorhombic 
phase with {a — b)/{ab) ^ 0.3%, large orbital-splitting 


\Exz — Eyz\ of order 50 meV is observed at X,Y-points by 
ARPES studies in BaFe 2 As 2 [30], NaFeAs [31], and FeSe 
[32-40]. Especially, noticeable deformation of the Eermi 
surfaces (ESs) with C 2 -symmetry is realized in EeSe, be¬ 
cause of the smallness of the Eermi momenta. In EeSe, 
Ref. [38] reports that the orbital splitting Exz — Eyz is 
positive at F-point, whereas it is negative at X,Y-points. 
This sign-reversing orbital splitting is not realized in the 
non-magnetic orthorhombic phase in NaFeAs [31]. In 
addition, the e-FSl at X-point is deformed to two Dirac 
cone Fermi pockets in thin-film FeSe [37, 40]. The aim of 
this study is to explain these nontrivial electronic states 
in the orbital-ordered states based on the realistic multi¬ 
orbital Hubbard model. 

Microscopically, the orbital order is expressed by the 
symmetry breaking in the self-energy. In the mean-field 
level approximations, however, the self-energy is constant 
in A;-space unless large inter-site Coulomb interactions 
are introduced [13]. For this reason, we have to study 
the non-local correlation effect beyond the mean-field 
theory, based on the realistic Hubbard model with on¬ 
site Coulomb interaction. We will show that the strong 
positive feedback between the nematic orbital order and 
(72-symmetric spin susceptibility plays the essential role. 

In this paper, we study the origin of the orbital or¬ 
der in Fe-based superconductors, by considering the co¬ 
operative symmetry breaking between the self-energy 
and spin susceptibility self-consistently. Experimen¬ 
tally observed strong /c-dependent orbital polarization 
{AExz{k), AEyz{k)) is given by the non-local self-energy 
with (72-symmetry. In the EeSe model, we obtain the 
sign-reversing orbital-splitting Exz—Eyz between F-point 
and X,Y-points reported in Refs. [38]. In addition, two 
Dirac-cone Fermi pockets emerge around X-point when 
the Coulomb interaction is larger [37, 40]. Thus, impor¬ 
tant key experimental electronic properties in the orbital- 
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ordered phase are satisfactorily explained. 

Hereafter, we denote the five d-orbital 
dyzi dxyi dx 2 _y 2 as / = 1, 2, 3,4, 5. We study the realistic 
eight-orbital d-p Hubbard models 

Huir) = i^M + (M = LaFeAsO and FeSe), (1) 

where is the first-principles screened Coulomb po¬ 
tential for d-orbitals in Ref. [41]: The averaged Coulomb 
interaction U = | = 7.21 eV for FeSe is much 

larger than U = 4.23 eV for LaFeAsO, since the number 
of the screening bands, by which U is reduced, is small 
in FeSe [41]. In contrast, the averaged Hand’s coupling 
J = ^ Sz>m similar in all compounds since the 

screening on the Hand’s coupling is small. For this rea¬ 
son, the ratio JjU = 0.0945 in FeSe is much smaller than 
the ratio J/U = 0.134 in LaFeAsO. The factor r(< 1) is 
introduced to adjust the spin fluctuation strength. 

in Eq. (1) is the first-principle tight-binding model 
introduced in Ref. [27]. The band dispersion is the 

solution of dei{z~^e^-ii — hPj^{k)) = 0, where is the 

kinetic term and is the diagonal mass-enhancement 
factor: = l/zi • Sip. For LaFeAsO, we put 

z~^ = 1. For FeSe, we put l/z^ = 1.6 and Xjzi — 1 
(/ 7 ^ 4) to represent the strong renormalization of the 
dxy-OThital band [32]. Figures 1 show the bandstruc- 
tures and the FSs in the FeSe and LaFeAsO models with 
ntot = 2 ^ fi^k) = 12. In FeSe, each Fermi pocket is 

very shallow [42]. The detail explanation for is given 
in the Supplemental Material (SM) A [43]. 



JL 

- ^ ^ 

e-FS2 

1 Ttj 

7 - Oj - 

^-FS2 


^-FSl X 


h-FSl 

ky 0 


kyO 

(g Q: 


h-FS2 e-FSl 


h-FS2 e-FSl 

-TT 

r\ 

-K 

, O 


kx kx 


FIG. 1: (color online) Bandstructures of the eight-orbital 
models for (a) FeSe and (b) LaFeAsO in the unfolded Bril- 
louin zone at T = 50 meV. FSs for the FeSe and LaFeAsO 
models are shown in (c) and (d), respectively. The colors 
correspond to 2 (green), 3 (red), and 4 (blue), respectively. 

In this paper, we calculate the /^-dependence of the 
self-energy using the self-consistent one-loop approxima¬ 


tion, which have been applied to various single-orbital 
models [45] and multi-orbital models [12] in literature. 
The Green function in the orbital basis is 

G{k) = {z~^i€n + fi- huik) - AE(fe))“\ (2) 

where k = (/c, = (2n + l)7rT) and AE(/c) is symme¬ 

try breaking self-energy, and z~^ is the diagonal mass- 
enhancement factor. The self-energy in the one-loop ap¬ 
proximation [12, 45] is given as 

E;,p(fc) = E[^;, + r y] Vi^m-,l',m'{q)Gm,m'{k - q), (3) 

q,m,m' 

where E^, = - is the Hartree 

term: Anm,m' = , 0 ) ~ ,(j)o- The 

non-local interaction V(q) in Eq. (3) is given as 

V(9)r* + (9)f“ -1 [f(g)f^ + f*x°(9)f* - 

+ T'^)], where = V(9)(l - 

f5,c^0(^))-l = -TT.kGl,m(k + 

q)Gm',i'{k). Here, is the bare Coulomb interaction 
for the charge (spin) channel given in the SM A [43]. 

Erom Eq. (3), the Hi^-type symmetry breaking self¬ 
energy included in Eq. (2), which is orbital-diagonal, is 
derived as 

AE,(fc) = Re IE,,((fc, e„) - Efe„)| (4) 

for ; = 2,3, where Ef^*(fc) = {Hi^iik) + E5-i,5-Kfc'))/2 
(k' = (ky^kx)) is the Aig-component of the self-energy. 
In the present study, we calculate Eqs. (3)-(4) self- 
consistently [46]. We will show that AE/ emerges due 
to the strong positive feedback between the nematic or¬ 
bital order and (72-symmetric spin susceptibility: That 
is, near the nematic transition, infinitesimally small ne¬ 
matic orbital order enhances the spin susceptibility at 
q = (7r,0), and the enhanced spin-fluctuation-mediated 
interaction V (g) in turn enlarges the orbital order. 

Eirst, we study the EeSe model, using 64 x 64 k- 
meshes and 512 Matsubara frequencies at T = 50 meV. 
In EeSe, the orbital order with AExz{k) = AT^ 2 {k) and 
AEyz{k) = AT>s{k) emerges when the spin Stoner fac¬ 
tor as, which is the maximum eigenvalue of r^x^(gr), 
is larger than 0.82 (r > 0.253). The magnetic or¬ 
der is realized when = 1. In Eigs. 2 (a) and 
(b), we show respectively the obtained orbital polar¬ 
ization {AExz{k), AEyz{k)) and the spin susceptibility 
X’iq) = 'Ei,mXli-,m,miq) for = 0.868. Due to the 
positive AEyzpC), the e-ESl around X-point is modi¬ 
fied to the pair of the Dirac-cone Eermi pockets. An¬ 
other electron-pocket around Y-point, e-ES2, is enlarged 
by the negative AExzC^)- We stress that AExz{yz) (fe) 
changes its sign along the (a,)-axis shown by the broken 
line in Eig. 2 (a). Due to this sign reversal, the outer 
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FIG. 2: (color online) (a) Orbital sign-reversing polarization 
AEyz{k)) and (b) x^{q) the FeSe model for 
as = 0.868 (r = 0.257). The realized FSs are shown by the 
green lines. The splitting Exz — Eyz at F-point is positive, as 
observed in Ref. [38]. (c) The polarization at F-, X- and Y- 
points. (d) An = rixz —riyz and (e) as for 0.261 > r > 0.251. 


hole-pocket (h-FS2) is elongated along the /c^-axis, as 
observed experimentally [38]. Due to the orbital polar¬ 
ization, as increases and x^{q) shows the C 2 anisotropy 
X^(7r,0) > X^(0,7r) shown in Fig. 2 (b) [6, 47]. 

In Fig. 2 (c), we show the r-dependences of the or¬ 
bital polarization at F-, X- and Y-points. With increas¬ 
ing r, the spin-fluctuation-driven orbital order appears 
as a second-order transition at = 0.82. The rela¬ 
tions AExz(y) < 0 and AEyz(X.) > 0 hold in the or¬ 
dered state. In FeSe, the orbital splitting Exz — Eyz 
at F-point is positive, so h-FS2 is elongated along the 
ky-diX.is. Such sign-reversing orbital order does not oc¬ 
cur in the LaFeAsO model (see Fig. 5). Since dxz- and 
dyz-oihitsih are exchanged by 7r/2-rotation, the obtained 
order {AExz{k)^ AEyz{k)) belongs to the Big representa¬ 
tion (=(i-wave), in spite of AExz{k) 7^ AEyz{k) [13, 36]. 

Figure 2 (d) shows the difference An = Uxz — '^yz- 
An ^ +10“^ will induce the small lattice deformation 
(a — 6)/(a + 6) ^ 0.2% in FeSe due to small e-ph in¬ 
teraction. When An % 0, the spin Stoner factor as is 
strongly enlarged as shown in Fig. 2 (e), consistently 
with the enhancement of spin fluctuations observed by 
NMR [20, 21] and neutron [22, 23] measurements. For 
a fixed r, An shows a mean-field-type second-order T- 
dependence below Tstr [27], 

In Fig. 3, we display the C 2 bandstructures obtained in 
the FeSe model for (a) as = 0.854 and (b) as = 0.868. In 
both cases, AEyz{kx^^) and AExz{^^ ky) show sign rever¬ 
sal, and Exz — Eyz > 0 at F-point. The shape of the FSs 
in Fig. 3 (c) for as = 0.854 (Fig. 3 (d) for as = 0.868) is 
consistent with the FSs observed in bulk FeSe (thin-fllm 
FeSe). This result indicates that the electron correlation 
in thin-fllm FeSe is slightly stronger, consistently with 



FIG. 3: (color online) Orbital polarizations and bandstruc¬ 
tures in the FeSe model along Y ^ F ^ X for (a) as = 0.854 
and (b) as — 0.868. The corresponding FSs at T = 50 meV 
are shown in (c) and (d), respectively. 


the higher Tgtr ~ 120K in thin-fllm FeSe. 

In real FeSe, the inner pocket sinks under the Fermi 
level due to the spin-orbit interaction. Similar single 
hole-FS model is introduced in the SM B [43] by shifting 
the dxy-oihiidl level at F-point. It is verified that the 
sign-reversing orbital polarization emerges in the FeSe 
model with single hole-FS in the SM B [43]. 

Here, we explain that the origin of the orbital order 
is the positive feedback between the orbital polarization 
{AExz{k), AEyz{k)) and the spin susceptibility. For this 
purpose, we introduce the following two simplified self¬ 
energies and perform the self-consistent calculations: 

SAL(fc) = SH,+T ^ Vi,ra;l',ra'{q)G^,^?{k-q)ib) 

q,m,m' 

q,m,m' 

where the superscript “AE = 0” means the absence of 
the C 2 orbital polarization. In Eq. (5) (Eq. (6)), 
only the feedback effect from the symmetry-breaking 
spin susceptibility (Green function) is included. 
contains the Aslamazov-Larkin term (Maki-Thompson 
term) that is the second-order (first-order) term with 
respect to in Eig. 4 (a), which shows the expan¬ 
sion of the self-energy in Eq. (3) by using the relation 
G = + G^=^ • AE • + 0(AE2); see the 

SM C [43] for detailed explanation. Eigure 4 (b) shows 
the orbital polarizations AEy^ and at X-point 

derived from Eqs. (5) and (6), respectively, as functions 
of r at T = 50 meV. Here, AE^^^ is the orbital polar¬ 
ization derived from Eq. (3), shown in Eig. 2 (c). The 
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closeness of and AEf}} means that the orbital or- 

der is mainly driven by the Aslamazov-Larkin term 
Figure 4 (c) shows that 0^5 is strongly enlarged by the 
orbital order due to similarly to as in Fig. 2 (e). 

In contrast, as in Fig. 4 (d) is suppressed because of the 
wrong q = (7r,0) nesting on the dyz-oihital due to 
Thus, the origin of the orbital order is the strong positive 
feedback between the orbital polarization and C 2 spin 
susceptibility described by the Aslamazov-Larkin term 
in Eq. (5). 

Although is the main driving force of the orbital 
order, the hole-pocket is elongated along the kx-axis in 
the ordered state driven by as shown in the left 

inset of Fig. 4 (b). Therefore, the Maki-Thompson term 
is indispensable to reproduce the sign-reversing orbital 
polarization; see the right inset of Fig. 4 (b). In SM 
C [43], we explain analytically why the Maki-Thompson 
term gives the sign reversal. Thus, although AE^^{kx,0) 
is always positive as shown in the left inset of Fig. 4 
(b), AEy^^{kx,0) in Fig. 2 (a) given by Eq. (3) changes 
to negative at ^ 0, due to the sign-reversing Maki- 
Thompson term that is small in magnitude at as ^ a^^. 



FIG. 4: (color online) (a) Expansion of the self-energy with re¬ 
spect to AE. The MT term (AL term) is included in |]MT(al) 

(b) r-dependences of AEy^ and AE^^ at X-point. AEy^^ is 
the orbital polarization derived from Eq. (3) shown in Fig. 2 

(c) . (c) as due to and (d) as due to 

Next, we show the numerical results for the LaFeAsO 
model with z~^ = 1. Then, the orbital order is realized 
for as > 0.91. In Fig. 5 (a), we show the obtained or¬ 
bital polarization at as = 0.98 and T = 47 meV. The 
e-FSl around X-point is smaller than e-FS 2 around Y- 
point due to the orbital polarization AExz{y) < 0 and 
AEyziA) > 0. In addition, AEyz{kx,t)) {AExz{^-,ky)) 
is always positive (negative), and therefore the outer 
hole-pocket is elongated along the kx-ax.is. This result 


is consistent with the Fermi surface deformation in the 
non-magnetic orthorhombic phase in NaFeAs [31]. In 
the ordered state, strong in-plane anisotropy of the spin 
susceptibility X^(7r,0) ^ X^(0,7r) is realized as shown in 
Fig. 5 (b) [47], consistently with the neutron scattering 
experiments. The bandstructure is shown in Fig. 5 (c). 
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FIG. 5: (color online) (a) Obtained sign-preserving orbital 
polarization [/\Exz{k)^ AEyz{k)) in the LaFeAsO model for 
as — 0.98 (r = 0.376). The FSs are shown by the green lines. 
The corresponding spin susceptibility and the bandstructure 
are shown in (b) and (c), respectively. 


In this study, we succeeded in explaining (i) the small 
critical Stoner factor a'g and (ii) the relation Exz > Eyz 
near F-point in FeSe. An origin of (i) is the smallness 
of the ratio J/U^ as verified by the SC-VC theory [27] 
and the renormalization group (RG) theory [48, 49]. An¬ 
other origin of (i) is the smallness of the FSs in FeSe, 
since the three-point vertex A 3 = Sx^{q)/SAE in the 
Aslamazov-Larkin term increases when the particle-hole 
asymmetry is large [27]. The relation (ii) can be realized 
by the sign-reversing Maki-Thomson term. Although the 
Maki-Thompson term is usually small, (ii) is actually re¬ 
alized when the FSs are smaller since the minimum of the 
Aslamazov-Larkin term AEy^{k) shifts to k = 0. Re¬ 
cently, the advantage of the small FSs for the nematicity 
had been stressed by the RG study in Ref. [50]. 

In summary, we investigated the C 2 symmetry break¬ 
ing in the self-energy and susceptibility self-consistently, 
and explained the experimental C 2 -symmetric FSs and 
X^{q)' III fhe FeSe model, experimental deformation 
of the FSs due to the sign-reversing orbital polarization 
is satisfactorily reproduced. In the LaFeAsO model, in 
contrast, spin fluctuations are strongly enlarged by the 
sign-preserving orbital polarization. Thus, the key ex¬ 
periments below Tstr are satisfactorily explained. The 
orbital order originates from the positive feedback be¬ 
tween the nematic orbital order and spin susceptibility 
due to the Aslamazov-Larkin term, and the sign reversal 
of AEyz{kx, 0) and AExz{ 0 , ky) in FeSe is caused by the 
Maki-Thompson term. 
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Sign-reversing orbital polarization in the nematic phase of FeSe 
due to the symmetry-breaking in self-energy 
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A: Eight-orbital models for FeSe and LaFeAsO, 
effect of the spin-orbit interaction 


We introduce the eight-orbital d-p models for FeSe and 
LaFeAsO analyzed in the main text. We first derived the 
first principles tight-binding models using the WIEN2k 
and WANNIER90 codes. Next, in order to obtain the ex¬ 
perimentally observed Eermi surfaces (ESs), we introduce 
the A;-dependent shifts for orbital /, SEi{k), by introduc¬ 
ing the intra-orbital hopping parameters, as we explain 
in Ref. [27]. We shift the dxy-oihitsl band [d^-^^/^^^-orbital 
band] at (F, M, X) points by (0, —0.25, +0.24) [(—0.24, 
0, +0.12)] for FeSe, and (0, 0, +0.05) [(-0.05, 0, +0.05)] 
for LaFeAsO. The unit is eV. 


Next, we discuss the bandstructure in the presence of 
the spin-orbit interaction (SOI), which is expressed as 
A li (Ti. The matrix elements of li are given in Ref. 
[44]. To study the effect of the SOI, we have to use the 
ten-orbital model since the “unfolding” is prohibited by 
the SOL In addition, the SU(2) symmetry for conduction 
electron spins is violated by the SOL For these reason, 
in the presence of the SOI, the numerical study of the 
self-energy calculation becomes very difficult. 
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FIG. SI: (color online) Bandstructures and hole-FSs in the 
FeSe model, (a) A = 0 (without SOI) and (b) A = 50 meV 
(with SOI) for Exz — Eyz = 0 at F-point. (c) A = 0 (without 
SOI) and (d) A = 50 meV (with SOI) for Exz — Eyz = 40 
meV at F-point. 


Figures SI (a) and (b) show the bandstructures and 
hole-FSs for A = 0 and A = 50 meV, respectively, in the 
FeSe model in the absence of orbital polarization. The 
main effect of the SOI on the hole-FSs is that the inner 
pocket shrinks or sinks below the Fermi level. In the 
presence of the orbital polarization Exz — Eyz = 40 meV 
at F-point, the bandstructures and hole-FSs are modified 
to Figs. SI (c) for A = 0 and (d) for A = 50 meV. In both 
cases, the outer hole-FS is mainly composed of the dxz- 
orbital, which is consistent with the experimental report 
[38], when we set the radius of the hole-FS just slightly 
smaller than the present model (/i = +10 meV). The 
hollows of the hole-FS near the /c^^-axis become smaller 
in the presence of the SOL 


channel in the main text is 
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Finally, the bare Coulomb interaction for the spin 


Also, the bare Coulomb interaction for the charge channel 
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Here, Ui^i^ U[ i, and J/,// are the first principles Coulomb 
interaction terms given in Ref. [41]. 
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FIG. S2: (color online) (a) (AExz{k), AEyz(k)) and (b) x^{q) 
in the FeSe model for as = 0.860 (r = 0.506) at T = 25 meV 
in the case of z = 1/2. The obtained C 2 -deformation of the 
FSs and that of x^{q) very similar to Figs. 2 (a) and (b) 
in the main text. 

In the main text, we put the renormalization factors 
as ( 2 ^ 1 , ^ 2 , ^ 3 , ^ 4 , 2 : 5 ) = (1,1,1,1/1.6,1). In the case of 
(^ 1 , 2 ^ 2 , 2 ^ 3 , 2 : 4 , 2 : 5 ) = 2 ;( 1 , 1 , 1 , 1 / 1 . 6 , 1 ), we proved that in 
Ref. [27] that essentially same numerical results are ob¬ 
tained by replacing r and T with rjz and Tz^ respec¬ 
tively. To verify this scaling relation, we perform the 
numerical study for 2 ; = 1 / 2 : Figure S 2 shows the ob¬ 
tained (a) {AExz{k), AEyz{k)) and (b) x^(gr) in the 
FeSe model for ^ = 1 / 2 , r = 0.506 and T = 25 meV 
{as = 0.860). These results are very similar to those for 
2 ; = 1 , r = 0.257 and T = 50 meV {as = 0.868) shown 
in Fig. 2 in the main text. Thus, the nematic orbital 
order is obtained by the present theory even if the band 
renormalization effect is taken into account. 

In the present study, we neglected the damping of the 
quasiparticle given by the imaginary part of the self¬ 
energy. We expect that this effect becomes small at suffi¬ 
ciently low temperatures, around the realistic Tgtr ^100 
K. 


B: Orbital polarization in the single-hole-pocket 
FeSe model 

In FeSe, the main effect of the SOI on the hole FSs is 
that the inner pocket sinks under the Fermi level. Sim¬ 
ilar bandstructure and the FS topology are realized by 
shifting the d^^y-orbital level at F-point by —0.6 eV, as 
shown in Fig. S3. Based on this “single-hole-pocket FeSe 
model”, we perform the numerical study of the orbital 
ordered state. 


FeSe (single-hole-pocket model) 




FIG. S3: (color online) (a) Bandstructure and (b) FSs in 
the single-hole-pocket FeSe model, by shifting the dx 2 /-orbital 
level at F-point by —0.6 eV. 


In the present model, the nematic orbital order is ob¬ 
tained for as > 0.846 (r > 0.253). In Fig. S4, we show 
(a) the orbital polarization {AExz{k)^ AEyz{k)) and (b) 
C 2 spin susceptibility X^( 7 r, 0 ) > X^( 05 ' 7 r) obtained for 
as = 0.877. Due to the orbital polarization, the size of 
the e-FSl around X-point is reduced, and two Dirac-cone 
Fermi pockets appear. Another electron-pocket around 
Y-point, e-FS 2 , is enlarged. In addition, the single hole- 
pocket is elongated along the /c^-axis. Such FS deforma¬ 
tion is essentially the same as that in the original FeSe 
model in the main text. In Figs. S4 (c)-(e), we show the 
r-dependences of the orbital orders. An = Uxz—nyz, and 
as^ respectively. 

In Fig. S5, we display the obtained bandstructures 
and orbital polarization for (a) as = 0.865 and (b) as = 
0.877. In both cases, AExz{y) — ^Eyz(X.) is negative, 
whereas AExz{^) — AEyz{T) is positive. The obtained 
FSs are also shown in Figs. S5 (c) and (d), respectively. 
These results are almost the same as the results in Fig. 
3 in the main text. 

Thus, the sign-reversing orbital polarization in FeSe, 
which is the main result of the main text, is also realized 
in the single-hole-pocket FeSe model. This result indi¬ 
cates that the sign-reversing orbital order survives in the 
presence of the SOI. 
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FeSe (single-hole-pocket model) 



FIG. S4: (color online) (a) Obtained orbital polarization 
{AExz{k), AEyz{k)) and (b) spin susceptibility in the single- 
hole-pocket FeSe model at as = 0.877 (r = 0.256). The real¬ 
ized FSs are shown by the green lines, (c) The orbital polar¬ 
ization, (d) An = Uxz—riyz, and (e) as for 0.258 > r > 0.251. 


FeSe (single-hole-pocket model) 



FIG. S5: (color online) Orbital polarizations and bandstruc- 
tures in the single-hole-pocket FeSe model along Y ^ F ^ X 
for (a) as = 0.865 and (b) as = 0.877. The corresponding 
FSs are shown in (c) and (d), respectively. 


C: Comment on the self-consistent vertex correction 
(SC-VC) method 

In Ref. [27], the present authors studied various model 
Hamiltonians for Fe-based superconductors using the SC- 
VC theory. In the FeSe model with =3, The critical 


spin Stoner factor for the orbital order is = 0.86 
at T = 50 meV. In the LaFeAsO model with = 1 
without the d-orbital level shifts, = 0.975. These 
critical spin Stoner factors are slightly larger than those 
obtained by the C 2 self-energy study in the main text. 
In FeSe, the orbital order is realized even when the spin 
fluctuations are weak in FeSe, due to the smallness of 
the ratio J/U and the broadness of the da-^^/^^^-orbital 
spin susceptibility, as explained in Ref. [27]. 

Here, we present an analytic explanation that the 
present self-consistent C 2 self-energy analysis is essen¬ 
tially equivalent to the SC-VC theory. In the main text. 
Fig. 4 (a) shows the diagrammatic expression of the self¬ 
energy in Eq. (3), by expanding the right-hand-side with 
respect to AEi = AE/: (i), (ii), and (iii), (iv) corre¬ 
sponds to the Hartree term, the Maki-Thompson (MT) 
term, and the AL-terms, respectively. In deriving Fig. 
4 (a), we used the relation G = + G^^^ • AE • 

gAE=o ^ 0{AE‘^). Based on Fig. 4 (a), the linearized 
self-consistent equation for AEi (/ = 2,3) is obtained by 
introducing the eigenvalue Aorb as 

XorbAEiik) [U{Gi{k)}^AEi{k) 

k 

+{2U'-J){G5-i{k)yAE5-i{k)] (S3) 

+TElVi%q){Gi{k + q)yAEi{k + q) 

q 

+Tj2{Gl{k + q) + Gi{k- q))^{Vf{q)}^Ai{AE-,q), 

q 

where the first two lines give the intra-orbital and inter¬ 
orbital Hartree terms, respectively, and the third and 
fourth lines give the MT and AL terms, respectively. 
Here, we assumed the orbital-diagonal G and for 
AE = 0, and Ai{AE-,q) = -Tj:^{Giip)}^Giip + 
q)AEi(jp) is the three-point vertex. The eigenvalue Aorb 
exceeds unity when the orbital order is realized. 

First, we explain that the orbital order mainly 
originates from the AL term: We put AExz{k) = 
—AEyz{k) = AE to simplify the discussion, and multiply 
both sides of Eq. (S3) by —T^f^{Gi{k)}‘^ by dropping 
the MT term. Then, we obtain 

AorbX°(0) « {2U' -J- U){x°{0)y + V^(0), (S4) 

where x^(0) is the bare-bubble and A^(0) is the charge- 
channel AL-VC for / = 2 or 3. Thus, the orbital or¬ 
der (Aorb = 1) is realized when A^(0) > x^(0) — {U — 
5J){x^( 0)}2 - 5J/7/2 for U = U' E 2J, if we put 
X^(0) 1/U for a qualitative estimation. This result 

is consistent with the condition for ac = 1 in the SC-VC 
theory A^(0) = (7/-5J)-^-x^(O) - bJ/U‘^ for J/U < 1 
[27]. Thus, the present theory is a natural extension of 
the SC-VC theory for the orbital-ordered state. 

We also explain that the MT term gives the sign- 
reversing orbital polarization in EeSe as explained in the 
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main text. Considering the fact that has the 

peak at (7r,0), the AE-linear MT term in Eq. (S3) for 
dyz-oihitsl at ~ 0 is approximately given as 

- E Vsf,^°iq)diEk+^ - i^)AEy,{k + q) 

Q 

^ - |E 0XS5) 

where is the ^z-orbital band along E-X shown in 
Fig. 1 (a) in the main text, and we used the relation 
^ Yjni^^n E 1^- - /i) at T ^ 0. There¬ 

fore, the MT term tends to induce the sign reversal of 


AEyz(xz){k) along the kx(y)-8iX.is, as depicted in the right 
inset of Fig. 4 (b) in the main text. Since AEy^ — AE^^ 
is the smallest at F-point in FeSe, the sign-reversing or¬ 
bital polarization is realized by the MT term. 

In the present C 2 self-energy study, the higher-order 
MT- and AL-terms are included, and therefore the ob¬ 
tained is slightly smaller than that given by the SC- 
VC theory. On the other hand, the feedback from the 
enhanced orbital susceptibility given by the VC is absent 
in the C 2 self-energy study. We will explain the relation 
between the C 2 self-energy theory and the SC-VC theory 
in more detail in future publications. 



